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ABSTRACT 

The expansion of ionization fronts in uniform and sphericaUy symmetric power-law density distri- 
butions is a well-studied topic. However, in many situations, such as a star formed at the edge of 
a molecular cloud core, an offset power-law density distribution would be more appropriate. In this 
paper a few of the main issues of the formation and expansion of HII regions in such media are outlined 
and results are presented for the particular cases where the underlying power laws are and . 
A simple criterion is developed for determining whether the initial photoionized region will be un- 
bounded, which depends on the power-law exponent and the ratio of the equivalent Stromgrcn radius 
produced by the star in a uniform medium to the stellar offset distance. In the expansion stage, the 
ionized volumes will eventually become unbounded unless pressure balance with the external medium 
is reached before the ionization front velocity becomes supersonic with respect to the ionized gas. 
Subject headings: H II regions — hydrodynamics — ISM: kinematics and dynamics 



1. INTRODUCTION 

The theoretical study of the formation and evolution of 
H II regions in non uniform media is motivated by the ob- 
servational result that the molecular clouds where mas- 
sive star form ge nerally possess density gradients (e.g., 
iLada et al]l2007| ). The assumed density law is taken to 
be a power law n oc r~" with exp onents ranging from 1 
to 3 ( Arquilla fc Go ldsmith' I T985f ) obtained from molec- 



ular line studies. Submillimeter continuum dust emis- 
sion observations of the molecular material around em- 
bedded ultracompact H II regions gives an exponent i n 
the range 1.25 to 2.25 (jHatchell fc van der Talj l2003f ). 
The most commonly found value for the power-law ex- 
ponent is a = 2, which correponds to an isothermal, 
self-gravitating sphere. However, steeper density laws, 
a < 4 have been inferred from radio conti nuum spectra 
of a sample of ultracompact H II regions (|Franco et al.l 
[2000[ ). 

Numerical studies of H II region expansion near the 
edge of a molecular cloud (the so-called "champagne 
flows"), were carried out by Tenorio-Tagle and coUab- 
orators (iTenorio-Ta gle 197a iBodenheimer et all 119791 : 
iBediin fc Tenorio-Taglelll98lh . In this model, an H II 
region formed near the edge of the cloud "breaks out" 
into the low density, intercloud medium during the ex- 
pansion stage and the ionized gas can achieve supersonic 
velocities of around 30 km s~^ in the acc elerating flow. 
This p roblem was studied in more detail bv lHennev et al.l 
(|2005l ). and applied to the Orion Nebula H II region. 
They found that a long-lived quasi-stationary phase ex- 
ists in which the ionized flow is approximately steady in 
the frame of reference of the ionization front, during this 
phase, the flow structure is determined entirely by the 
distance of the ionizing star from the front and the cur- 
vature of the front. The curvature determines whether 
the flows are more "champagne like" or "globule like". 
The curvature of the ionization front depends on the lat- 
eral density distribution in the neutral gas. 

^ Email: j.arthur@astrosmo.unam.mx 



[Franco et al.l ()1989l I1990D developed analytical solu- 
tions for the formation and expansion of H II regions in 
spherically symmetric, power-law density distributions, 
where the ionizing star is assumed to be at the center 
of the density distribution. The main result of their pa- 
per is that for spherical clouds with power-law density 
distributions r"", there is a critical value of the expo- 
nent, ttcrit = 3/2 above which the cloud becomes com- 
pletely ionized during the expansion phase of the H II 
region. For power-law exponents higher than this crit- 
ical value, two regimes were identified: a slow regime, 
in which the expansion of the dense cloud core is mildly 
supersonic and has almost constant velocity, correspond- 
ing to 3/2 < a < 3, and the fast regime, where the core 
expansion drives a strongly accelerating shock, which 
correspon d s to a > 3. This problem was revisited by 
IShu et aTl (j2002f ). who developed a self-similar solution 
for the internal motions in the ionized gas in the a > 3/2 
c ases. 

iFranco et all (|1989L Il990( l also considered self- 
gravitating isothermal disk density distributions where 
the density fall-off is steeper than r~^^^. In this case, 
the H II region becomes unbounded in a conical section 
of the disk, where the opening angle of the cone depends 
on the ratio of the Stromgren radius in the midplane den- 
sity to the density distribution scale height. During the 
expansion stage of such H II regions, the flattening of the 
internal ionized gas density distribution can lead to the 
trapping of an initially unbounded ionization front. 

More recently, photoionized regions in a variety of ax- 
isymmet ric density d istributions were modeled numeri- 
cally by [Arthur fc Hoard (j2006l ) , who also included the 
central star's stellar wind. In that work, it was found that 
the hot bubble created by the shocked stellar wind results 
in a champagne flow being set up outside the swept-up 
stellar wind shell, leading t o mixed kinemat i cs fro m an 
observational point of view. [Arthur fc Hoard (|2006D also 
considered the possibility of stellar motion up the den- 
sity gradient, which introduces the added complication 
of bowshock kinematics to the problem. 
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iMac Low et al.l (|2006[ ) performed low-resolution three- 
dimensional numerical simulations of H II region dynam- 
ical evolution in a collapsing molecular cloud. Here, the 
density distribution had a core-envelope structure, where 
the envelope had a power-law (r~^) density distribution. 
Both turbulent and non-turbulent cases were modeled. 
In the non-turbulent case, even for ionizing sources off- 
center from the core, the H II region expansion was nearly 
spherical for the parameters adopted and the timescale 
of the simulation. 

The photoionized region s produced by mov ing stars 
were studied numerically bv lFranco et al.l (|2007f ) . In that 
paper, the H II regions produced by a star off-center in 
a spherically symmetric density distribution were mod- 
eled, where a range of stellar velocities (between and 
12 km s"'^) were considered. It was found that the H II 
region produced by a star off-center in an underlying 
density distribution did not become unbounded, while 
that in a r^^ density distribution did, for the case of a 
stationar y star. This appe ars to contradict the main re- 
sult of Franco et al.l (|1990( ). This result is the motivation 
for the present work, where we aim to provide simple 
tools for analyzing the formation and expansion of H II 
regions off-center in spherically symmetric density distri- 
butions. 

The structure of this paper is as follows: in §[5] we out- 
line the problem and discuss the initial formation stage of 
an H II region off-center in a spherically symmetric den- 
sity distribution and find a criterion for the H II region 
to remain bounded. In §[3] we investigate the expansion 
stage and derive a simple differential equation that de- 
scribes the expansion of the ionization front along the 
symmetry axis. We assess the validity of our assump- 
tions about the internal structure of the ionized gas by 
comparing our analytical model with the results of nu- 
merical simulations in § |4l Finally, in § [5] we summarize 
our results. 

2. DESCRIPTION OF THE PROBLEM 

The underlying, spherically symmetric density distri- 
bution is described by 



n(r) 



(1) 



where Uc is the number density at the position of the 
star, Tc. We now change our reference system and rewrite 
the density distribution in axisymmetric coordinates cen- 
tered on the star position, with the distance from the 
star, R, found from 



- 2rci?cos(7r 



, 



(2) 



where 9 is measured from the symmetry axis and the 
center of the density distribution is in the direction 9 = tt, 
see Figure [1] The density distribution of Equation [U in 
terms of this new reference system, (i?, 9), becomes 



n(i?, 9) — ric 



2rrRcos9 



-q/2 



(3) 



Along the symmetry axis, 9 = 0, this can be written 

niy)=nAl + yr" , (4) 
where we have defined y = R/tc- 




Fig. 1. — Geometrical configuration. The center of the (spheri- 
cally symmetric) power-law density distribution is denoted by the 
open circle, the star is represented by the filled circle, a distance 
Tc from the density center. For a parcel of gas at X, distance r 
from the density center, the star-centered coordinates are {R,6), 
where R is the distance from the star and 8 is measured clockwise 
from the extension of the line joining the star to the density cen- 
ter. Rotational symmetry is assumed around the axis joining the 
star to the center of the density distribution. The position of the 
ionization front in this coordinate system is Ri. An approximate 
shjipe for the ionization front is also shown. 

For y < 1, expansion to linear terms in a Taylor series 
shows that 

n-iy) ~ nc{l - ay) , (5) 

which is no longer a power law. For y > 1 we find 
« Ticy^"", that is, the offset density distribution re- 
sembles the original power-law density distribution. Con- 
sequently, if the initial Stromgren radius, Rs, is smaller 
than the stellar offset distance, then we expect the ini- 
tial H II region to remain confined, even for power laws 
steeper than a > 3/2. 

2.1. Condition for the initial Stromgren region 

We can develop a quantitative criterion for decid- 
ing whether the initial photoionized region will remain 
bounded in an offset power-law density distribution by 
balancing recombinations and ionizations, assuming pure 
hydrogen, along the symmetry axis 9 = 0: 



y 



(6) 



where as is the case B recombination coefficient. This 
equation can be written 



(l-hy')"'"y" dy 



1 



(7) 
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and momentum conservation 



Fig. 2. — Ratio of ionization front radius along 9 = to uniform 
density Stromgren radius, y/ysc, as a function of uniform density 
Stromgren radius, j/sc, for power-law density distributions with ex- 
ponent cx = 2 (solid line) and a = 3 (dashed line). The asymptotes 
are i/ac = 1 for a = 2 and i/bc = 0.464 for a = 3. 

wfiere j/sc is tlie non-dimensionalized Stromgren radius 
for a uniform medium of number density ric, i.e., 



Vsc = 



1/3 



(8) 



In the "uphill" direction, = tt, the equivalent non- 
dimensionalized ionization balance equation is 
ry 



(1-2/') 



-2a 12 1 I 

V dy 



1 



(9) 



Equations [8] and [9] define integral equations for the ex- 
tent of the ionized region, y, along 9 = Q and 6* = tt, 
respectively, as functions of ysc and a. 
Integrating Equation [7] for a > 3/2 we obtain 



(1-2q)(2-2q)(3-2q) 
2 



-jy(l-[l+t/p-2") 



(10) 



In general, this yields a positive, real root for y (the non- 
dimensionalized ionization front radius) as long as 



1 3 2 

S^'^'^ (2a- l)(2a-2)(2a-3) 



(11) 



For example, when a — 2, the requirement is < 1, 
while for a = 3, the condition is y^^ < 0.1, i.e., ysc < 
0.464, in order that there be a real, positive root for y. 
In Figure [2] we plot the ratio y/ysc against for ysc 
values between 0.01 and 1.0 for the two examples above, 
a = 2 and a = 3. We see that for ysc less than the 
critical value, the value of the root y is not too different 
from ysc, but close to the critical value y/ysc — *■ oo- 

3. EXPANSION STAGE 

In order to study the expansion of the ionized gas after 
the initial formation of the photoion ized region we f ollow 
the same analysis as 'S pitzeii ()1968[) and iShul ()1992[ ). As 
the H II region expands, a neutral shock is sent out ahead 
of the ionization front into the ambient medium. 

Across the isothermal shock, the relevant equations are 
mass conservation 



PlsVls = P2sV2s , 



(12) 



Pis 



P2s + P2sV2s 



(13) 



where the subscript Is refers to preshock conditions and 
2s refers to postshock conditions and all velocities are in 
the frame of the shock. We also have the equation for 
the density jump across a weak-D ionization front 



P2i 

Pii 



1 

24 



1/2 



(14) 

where the subscript li refers to conditions conditions in 
front of the ionization front and 2i refers to conditions 
in the photoionized gas just behind the ionization front. 
Here, ci is the sound speed in the neutral gas and C2 is 
the sound speed in the ionized gas, and vu refers to the 
velocity of the gas ahead of the ionization front in the 
frame of the ionization front. 

Furthermore, we assume that the region between the 
shock front and the ionization front is thin, i.e. uu = U2s, 
where uu and U2s are the velocities of the gas ahead of 
the ionization front and behind the shock front in the 
frame of the star, and that the density of the ionized 
gas inside the photoionized region is spatially uniform 
but decreases with time, i.e. p2i = P2i(t). This is valid 
if the sound-crossing time of the ionized region is short 
compared to the expansion timescale. 

As a first approximation to the non-spherical shape of 
the ionization front in a non-uniform density distribu- 
tion, we assume that the radius of the ionization front 
can be described by 



Ri = clq + ai cos 9 — ao 
to first order in fi — cos 9, where 

R- + R+ 



aiP 



and 



ao = 



2 

R— — R-i 



(15) 



(16) 



(17) 



with i?_ corresponding to the shell radius when fj, = —1 
and R+ to the case p — I (see, e.g.. lDvsonlll977[ ). 
The volume contained within the ionization front Ri is 



V, 



27r 

^ J-i 

which in this approximation becomes 



al 



l-f 



(18) 



(19) 



Assuming uniform expansion, then the mass contained 
within a volume, V , interior to an internal shell of radius 
R must be conserved during the expansion. That is, 

P2iV = constant , (20) 

and the total volume of ionized gas satisifies the ioniza- 
tion balance equation, where the ionizing photon output 
of the star is assumed to be constant, that is 

p\iVi = constant . (21) 

Taking the time derivatives of Equations [20] and [21] we 
find that 

(22) 



dV 



2 dt 



4 



Arthur 



where we have evaluated the derivative for a shell located 
just inside the ionization front radius. Writing 



100 



dV 



dV dR 
'dR~dt 



and 



dV^ 
dt 



dV, dR, 
dRj dt 



(23) 



then we find that the gas velocity behind the ionization 
front is half the ionization front velocity, U2i — O.bUi, 
where U2i = ^1^. and Ui — This assumes that the 
expansion is uniform and that there are no lateral gas 
motions within the photoionized region. 

With the additional assumption that vu/ci <C 1, then 
equations [121 [131 [SI together with U2i = 0.5t7i, can be 
combined to give 



P2i 
Pis 



4c2 + Uf 



(24) 



where we have used the thin shell approximati on Uj Ki 
Us- T hese assumptions are the same as those of ISpitzeil 

dSH). 

The density ratio can be found from ionization balance 



(25) 



where Rq is the initial Stromgren radius formed in a 
medium of uniform mass density pc- We also have the 
offset power-law initial density distribution, written here 
for /i ± 1 : 



Pis = Pc 



(26) 



where is the offset distance of the star from the center 
of the density distribution and is the density at this 
point. Thus 



P2i 
Pis 



1/2 



rr 



in the directions /i = ±1 along the symmetry axis. 
Equation 1241 can be written in the form 



P2i 
Pis 



= 4 



P2t 
Pis 



(27) 



(28) 



where M — Ui/ 02- This constitutes a differential equa- 
tion for the radius of the ionization front, Ri, as a func- 
tion of time and can be solved num erically with standard 
techniques fe.g.. [Press et al.lll992f ). by considering direc- 
tions p = I and /i = — 1 together and using equations II 51 
[I6l[17l[l8]and[23 with initial conditions M = 2/^3 when 
f = 0, Ri = Ro and p2i = pis- Here, it is assumed that 
i?o is much less than the critical value found from solv- 
ing Equation 111! above, and so the initial ionized region 
radius {Ri at t = 0) is taken to be that for a uniform 
medium. 

In Figures [3] and [4] we can see the result of solving 
Equation [28] numerically for the cases a — 0, 2 and 3. 
In these figures, all distances are scaled to Tc and times 
are scaled with the sound crossing time of this distance 
in the ionized gas, r = rc/c2- The results presented 
are particular to the case Usc = Ro/fc = 0.01324, which 
corresponds to an ionizing source with Qh = lO''^ pho- 
tons s~^ and a density Uc = lO^/e cm~^, where we have 
also assumed that the case B recombination coefficient 
takes the value as = 2.6 x 10~^^ c m'^ s^^. These value s 
of the parameters are motivated by [Franco et al.l (|2007[ ) . 
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Fig. 3. — Radius of the ionization front (in units of the stellar 
offset radius, Vc) against time (in units of the sound crossing time, 
T = rc/c2). Top: fi = 1 direction. Bottom: fi = —1 direction. 
Density power laws are « = (solid line), a = 3/2 {dashed line), 
a = 2 {dot-dashed line), and a = 3 {dotted line). 

In the "upward" direction [p = —1, density increas- 
ing), as expected the radius grows more slowly in the 
a — 2 and 3 cases than for the uniform medium, clas- 
sical case a = 0, as shown in Figures [3] and [4l In the 
"downward" direction, p = 1, Figures [3] and [H clearly 
show that in the a = 3 case, the velocity starts to in- 
crease after about 5 sound crossing times and by 100 
sound crossing times the velocity of the ionization front 
has become supersonic with respect to the ionized gas. 
This means that the ionization front has overtaken the 
shock wave and will ionize out to infinity. In the a — 2 
case, although the velocity does start to increase again 
after about 13 sound crossing times, the increase is only 
very gradual, and within the timescale studied the ion- 
ization front does not become supersonic with respect to 
the ionized gas. 

In Figure[5]we show the ratio of the pressure in the ion- 
ized gas to that in the ambient medium, assuming that 
the latter is isothermal with a temperature of 100 K. 
From these figures we see that in the p — —1 direction 
the pressures equalize in about 1 sound crossing time 
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Fig. 4. — Velocity of the ionization front (in units of the ionized 
gas sound speed) against time (in units of the sound crossing time). 
Top: (1 = 1 direction. Bottom: /i = ~l direction. Linetypes same 
as in Fig. |3] 



for both a — 2 and a = 3. In the fj, — 1 direction, 
the a = 3 case shows that the pressures do not equahze 
before the ionization front begins to accelerate. In the 
a — 2 case, the pressures do equahze at around 7 sound 
crossing times for the adopted parameters (i.e., initial 
Stromgren radius, Rq). We remark that once the pres- 
sures become similar, the analysis of this section does not 
really hold, since the external pressure should be taken 
into account. However, we also comment that the fact 
that the pressures do become equal before the shock be- 
gins to accelerate in the a = 2 case means that, in prac- 
tice, the H II region will remain confined for this choice 
of initial Stromgren radius. 

In Figure [H] we show how the breakout time of the 
ionization front in the fJ. = 1 direction varies with initial 
Stromgren radius, for values of Rq/vc between 0.01 and 
1.0 for the case a — 2. As Rq/tc ^ I the position of 
the velocity minimum moves to earlier times and higher 
velocities. Also in this figure, we show how the internal 
to external pressure ratio varies with initial Stromgren 
radius. Only for Ro/rc S 0-02 is there a regime for which 
pressure balance could possibly halt the expansion and 
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Fig. 5. — Internal to external pressure ratio against time (in 
units of the sound crossing time). Top: fJ. = 1 direction. Bottom: 
fi = —I direction. Linetypes same as in Fig. [3] The horizontal line 
marks equal pressures. 

subsequent breakout of the H II region. 

4. NUMERICAL SIMULATION 

In order to test our assumption that the density within 
the ionized gas is approximately constant during the ex- 
pansion stage, we perform numerical simulations of the 
expansion of an off-center H II region in a power-law den- 
sity distrib ution. The density law for this simulation is 
taken from [Franco et all (|2007f ). namely 



n(r) = 



no exp 



noe 



-q/2 



for r < Vc, 
for r > Vr, 



(29) 



where no is the central density and Tc is the stellar offset 
distance from the center of the density distribution. At 
Tc the density has the value ric = n{rc) = no/e"^^. 

This density distribution is flatter for radii r < Vc than 
the simple power-law density distribution discussed in 
§§ [2] and [3] above. It is obtained by assuming that the 
acceleration due to gravity of the hydrostatic equilbrium 
solution in the halo, r > Tc, matches that of the core 
(defined as r < Tc) at r = Tc, with the added requirement 
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Fig. 6. — Top: Velocity of the ionization front (in units of the 
ionized gas sound speed) against time (in units of the sound cross- 
ing time) for different values of the initial Stromgren radius. From 
lowest line to uppermost line: Ro/tc = 0.01, 0.02, 0.05, 0.1, 0.2, 
0.5. Bottom: Ratio of internal to external pressures for the same 
initial Stromgren radius values. The horizontal line denotes equal 
pressures. 

that the acceleration due to gravity must be zero at the 
center of symmetry. This particular density law assumes 
that the acceleration due to gravity is linear with radius 
for r < Tc- We adopt this density distribution rather 
than a simple power law because the density does not go 
to infinity at the center of symmetry and the "core-halo" 
structure is more representative of real molecular clouds. 

Although we do not expect the ionization front ex- 
pansion of our numerical simulations to agree with the 
analytic expansion found in §[3] in the region between the 
stellar position and the center of the density distribution 
(i.e., r < Tc), because the density laws are very different 
here, we do expect to be able to make a direct compari- 
son in the halo region, r > rc, where the density laws are 
the same. Since we are primarily interested in whether 
the H II region remains bounded in the fJ- — 1 direction 
during the expansion stage, then it is sufficient that the 
density distributions are the same in this direction for 
radii r > Vc- 

The radiation-hydrodynamics code employed in 
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Fig. 7. — Logarithm of ionized gas density. Top: o = 2 numerical 
solution after 60,000 yrs. Bottom: a = 3 numerical solution after 
50,000 yrs. The black triangle marks the position of the star and 
the position (0, 0) denotes the center of the density distribution 
in each case. The heavy black curve represents the shape of the 
ionization front obtained from the approximation of Equation 1151 

these simulations has been described previously in 
[Arthur fc Hoard (|2006D . The model parameters specific 
to this work are a central density uq — 3.5 x 10^ cm^"^, 
and distance from density center to the star of Vc — 
0.15 pc. We assume isothermal, hydrostatic equilibrium 
in the initially neutr al 100 K gas, and include an exter- 
nal gravity field (see [Franco et al.ll2007i for details). The 
gravitational acceleration does not have any noticeable 
dynamical effect on the expansion of the H II region. 
The ionizing star produces 10"'* photons s~^. The grid 
size for these calculations is 250 x 600 cells, represent- 
ing a spatial size of 0.25 x 0.6 pc. For these parame- 
ters. Use = Rq/tc = 0.08235 for the a = 2 case and 
ysc — Ro/rc ~ 0.11492 for the a = 3 case. 

In Figures[71E]and[n]we show the results of calculations 
of off-center H II regions in a = 2 and a = 3 power-law 
density distributions. Figure[7]shows the ionized gas den- 
sity, i.e. the shape of the H II regions, for both power 
laws after 60,000 and 50,000 yrs respectively. The den- 
sity within the ionized regions is approximately constant. 
We note that for both models, the ionized volumes are 
approximately spherical, although the star is o ff-center. 
This was also noted by iMac Low et al.l (|2006f ) in their 
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Fig. 8. — Cuts along the fj, = ±1 axis of the a = 2 numerical 
simulation shown in Fig. [7] Top: density. Middle: velocity. Bot- 
tom: pressure. The thick solid line represents the fi = I direction 
and the thin solid line represents the fj. = —1 direction. Units are 
as indicated on the axes. 

three-dimensional numerical simulations. Overplotted 
on Figure [7] is the shape the ionization front would have 
according to the approximation of Equation 1151 We can 
see that the approximation slightly underestimates the 
volume of the ionized gas with respect to the numerical 
simulations. 

In Figures [8] and [9] we show cuts along the symmetry 
axis, i.e. the fj, = ±1 directions, for both cases. From 
these figures we see that the assumption of constant den- 
sity is reasonable for the ^ — 1 { "downhill" ) direction but 
not so good in the /i = — 1 direction, since the density in- 
creases with distance from the star. The pressure follows 
the same pattern as the density, as is to be expected in 
for an isothermal, T = 10* K H II region. Outside the 
ionized region, in both directions, the pressure is high be- 
cause the neutral shock sent out ahead of the ionization 
front pressurizes the gas. The velocities in the ionized 
gas are subsonic in the fi = —I direction and become 
supersonic in the /i = 1 direction. 

The ripples in the ii — —1 subsonic gas are due to 
the ionization front not being well resolved in this case, 
since the high densities lead to a very thin transition 
zone. The spike in the pressure is due to the heating 
and cooling being out of equilibrium at the unresolved 
ionization front. These numerical problems are not dy- 
namically important for these simulations, since we are 
interested in the expansion in the n — 1 direction, where 
the gas is supersonic and hence not affected by fluctua- 
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Distance from star [pc] 

Fig. 9. — Cuts along the = ±1 axis of the a = 3 numerical 
simulation shown in Fig. [7] Top: density. Middle: velocity. Bot- 
tom: pressure. The thick solid line represents the fi = 1 direction 
and the thin solid line represents the fi = —1 direction. Units are 
as indicated on the axes. 

tio ns in the thermal pr essure. This problem was noted 
bv iHennev et al.l (|2005f ). Their solution to the problem 
was to reduce the effective ionization cross section, (Tq, 
by a factor of 60. In this paper, we reduce ctq by a factor 
of 10. Reducing the ionization cross section makes the 
ionization front broader and so allows it to be resolved 
by the numerical scheme. Since the width of the ioniza- 
tion front is of the order of a few photon mean free paths, 
where Ip cx {nao)~^, then the ionization front is narrower 
in high density gas than in low density gas, and hence 
more of a problem to resolve numerically in the former. 
In the /i = 1 direction, the much lower densities mean 
that the ionization front is easily resolved by the numer- 
ical scheme. 

In Figure [Tol a comparison of the numerical simulation 
and the analytic solution for the same parameters in the 

= 1 direction is plotted. Here, the distance scale is in 
units of the Stromgren radius in a medium of uniform 
density ricC""/^, and time is given in units of the sound 
crossing time of the Stromgren sphere. Since the density 
distribution in the fj, = —1 direction of the numerical 
simulation is different to the analytical model discussed 
in Section [31 we do not show a comparison graph for 
this case, as it is meaningless. Instead, we focus on the 
fJ. = 1 direction, for which the numerical simulation and 
analytical model density distributions arc the same. We 
plot both the shock wave and the ionization front tra- 
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Fig. 10. — Radius (in units of stellar offset radius) against time 
(in units of sound crossing time) in the /i = 1 direction, showing 
ionization front radius {thick solid line), shock front radius {thin 
solid line), analytical solution for the same parameters {dot-dashed 
line) and uniform density analytical solution {dashed line). Top: 
a = 2 numerical simulation and corresponding analytical solution. 
Bottom: a = 3 numerical simulation and corresponding analytical 
solution. 

jectories, together with the analytical result. For times 
less than about half a sound crossing time, there is not 
a good agreement between analytical and numerical re- 
sults. This can be partly attributed to resolution effects, 
since the initial Stromgren radius in this high density 
medium is very small, and the details of the hydrody- 
namics are not particularly well resolved at this scale. 



For times greater than half a sound crossing time, the 
agreement between numerical and analytical solutions is 
very close for the remainder of the calculation. The ana- 
lytical solution falls between the shock front and ioniza- 
tion front paths for the whole of this time. 

We do not follow the numerical simulation beyond 
around 4 sound crossing times because part of the ionized 
gas volume leaves the computational grid at this point. 

5. SUMMARY 

In this paper we have developed an analytical treat- 
ment for the formation and expansion of an H II region 
off center in a steep power-law density distribution. We 
find that during the initial formation stage, the H II 
region will remain bounded as long as the ratio of the 
Stromgren radius in the equivalent uniform medium at 
the stellar position to the stellar offset position is less 
than a value that depends only on the power-law expo- 
nent (see Eq. llip . In the expansion phase, the ionization 
front will break out (turn supersonic with respect to the 
ionized gas) unless pressure balance between the internal 
photoionized gas and the external medium is achieved 
before this point. For Stromgren radius to stellar off- 
set radius ratios of more than 0.02, pressure balance is 
not possible for an isothermal 100 K medium in hydro- 
static equilibrium and so the H II region will become 
unbounded, eventually, in these cases. 

We performed numerical simulations to test the as- 
sumption of constant density within the ionized gas dur- 
ing the expansion phase adopted in § [3l We found that 
this assumption is valid in the "downhill" (6 — 0) di- 
rection from the star, but is not so good in the oppo- 
site ("uphill", 9 = it) direction. In the 9 — direction, 
the numerical solution and analytical solution follow the 
same expansion law once the elapsed time is greater than 
about half a sound crossing time. We also assessed the 
validity of the simple analytical approximation to the 
ionization front shape. Although the volume of ionized 
gas is slightly underestimated in the analytical approx- 
imation, the effect is not significant, as evinced by the 
agreement between the numerical and analytic expansion 
laws. This gives us confidence that our simple descrip- 
tion of the H II region expansion can provide insight into 
this complicated problem. 
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